Acoustic modes in metallic nanoparticles: atomistic versus elasticity modeling 



o 
o 

(N 



CZ3 



Nicolas Combf] 

CTVES ; CEMES (Centre d' Elaboration des Materiaux et d'Etudes Structurales); BP 94347, 
29 rue J. Marvig, F-31055 Toulouse, France, and 
Universite de Toulouse ; UPS ; F-31055 Toulouse, France 

Lucien SaviotQ 

Institut Carnot de Bourgogne, Universite de Bourgogne, UMR 5209 CNRS, 
9 Avenue A. Savary, BP 47870, F-21078 Dijon Cedex, France 
(Dated: July 6, 2009) 

The validity of the linear elasticity theory is examined at the nanometer scale by investigating 
the vibrational properties of silver and gold nanoparticles whose diameters range from about 1.5 
to 4 nm. Comparing the vibration modes calculated by elasticity theory and atomistic simulation 
based on the Embedded Atom Method, we first show that the anisotropy of the stiffness tensor in 
elastic calculation is essential to ensure a good agreement between elastic and atomistic models. 
Second, we illustrate the reduction of the number of vibration modes due to the diminution of the 
number of atoms when reducing the nanoparticles size. Finally, we exhibit a breakdown of the 
frequency-spectra scaling of the vibration modes and attribute it to surface effects. Some critical 
sizes under which such effects are expected, depending on the material and the considered vibration 
modes are given. 

PACS numbers: 63.22.-m, 63.22.Kn,81.05Bx 
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I. INTRODUCTION 

Describing structural and elastic properties at the 
nanometer scale has recently raised the interest of many 
scientists. Especially, the applicability of the linear elas- 
tic theory is questioned at these scales. Various ap- 
proaches can be used to test the validity of linear elas- 
ticity. Different groups have tried to measure the elastic 
properties, most often the Young modulus, by mechanical 
deformation while decreasing the system sizep],0,0,0,0]. 
Another approach to probe the elastic properties uses 
acoustic vibration of nanoparticles (NP) and is mainly 
based on light or neutron scattering. Acoustic vibration 
modes frequencies of NP measured by Raman scattering 
or pump-probe experiments are fairly well reproduced by 
linear elastic theory calculations even for NPs with sizes 
of a few nanometers d, 0, 0, 0] . 

The present work is concerned with the validity of 
elasticity within the vibration properties computed for 
systems whose sizes are smaller than the ones currently 
and experimentally explored. To achieve this goal, we 
compare the vibration modes calculated by linear elas- 
ticity and atomistic semi-empirical potential calculations 
in the case of metallic NP of diameter ranging from 1.4 to 
about 4 nm. The choice of metallic NP is motivated by 
the interest they arouse among the physicists community 
due to their plasmon properties opening many promising 
applications [l(| EH • 

Vibration properties of metallic NPs have been the 



'Electronic address: combe@cemes.fr 
^Electronic address: lucien.saviot@u-bourgogne.fr 



subject of many theoretical studies using the linear elas- 
ticity or atomistic simulation, and experiments. 

Since the nineteenth century, the linear elasticity has 
been used to calculate the vibrational modes of spher- 
ical particles [T^: such calculations generally agree well 
with experimental measurements of NP vibrations [71. 1 1 3j| . 
Different refinements of this model have been performed 
taking into account the anisotropy of the stiffness ten- 
sor, the non-spherical shapes of the nanoparticles and 
matrix effects [lj, [H, (3- The improvements of the ex- 
perimental energy resolution and of the synthesis of such 
NPs have allowed to demonstrate subtle effects due to 
the anisotropy of the stiffness tensor [ijj and due to the 
shape of the NPs 0,011. 

On the other side, using atomistic simulations, Raman 
and Kara[H, [2(| studied the vibrational density of state 
of metallic NPs. They showed that there is an enhance- 
ment in the vibrational density of states at low frequen- 
cies and an overall shift of the high frequency band be- 
yond the top of the bulk phonons when decreasing the 
NP size. The effect of the capillary pressure induced 
by the surface of the NPs has also been exhibited on 
the vibrational density of state using atomistic numeri- 
cal approaches [U [13]. Besides the study of the vibra- 
tional density of states, severals studies have focused 
on few given vibration eigenmodes of NPs. Focusing 
on Raman active vibration modes, the elasticity theory 
showed a very good agreement with atomistic calcula- 
tions in spherical germanium nanoparticle [23l l24j . Re- 
cently, a breakdown of frequency-spectra scaling of re- 
spectively, silicon and zinc oxide nanoparticles have been 
evidenced using atomistic simulations and attributed to 
surface effects 0, HI, H|. 

In this study, we would like to address three issues 
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concerning the vibration properties of metallic NPs com- 
paring the linear elasticity predictions and atomistic cal- 
culations. 

First, we compare the vibration modes of metallic 
NPs calculated using either the elasticity theory or an 
atomistic approach based on an Embedded Atom Model 
(EAM). A similar study has been performed by Cheng 
et al.[23| in the case of germanium NPs. However, 
their elastic calculations did not take into account the 
anisotropy of the stiffness tensor. In this study, we show 
that taking into account this anisotropy in the elastic 
calculations improves the agreement between both ap- 
proaches and that linear elasticity reproduces very well 
the vibration modes of NPs having a diameter of a few 
nanometers. 

Second, it is well known, that the vibration proper- 
ties of a NP containing N atoms can be described by 
3iV normal vibration modes. Therefore reducing the size 
of the NP also reduces the number of normal vibration 
modes. Elasticity, as for it, can not take into account this 
reduction because of its description using a continuum 
medium. Some studies include this reduction in elast icity 
theory by artificially limiting the number of modes [2 71]. 
In this study, we illustrate the reduction of the number 
of normal modes by showing the disappearance of certain 
modes while decreasing the NPs size. 

Third and finally, while elasticity predicts that the fre- 
quency of a given normal mode scales as the inverse of a 
characteristic length of the NP, surface effects can alter 
this law [3, HH]. We analyze this scaling in the case of 
silver and gold spherical NPs and exhibit a characteristic 
size under which the frequency scaling break down. The 
breakdown of this scaling law is related to the surface 
stress and surface relaxation effects. 

Section HT1 describes the numerical method used for the 
elastic calculations, the atomistic model and the proce- 
dure used for comparing the two models. Section IIIII 
addresses the three problems mentioned above. Most of 
the results presented here concern silver NPs, except for 
the last issue for which the case of gold NP has been 
considered. 



II. METHOD 



where p is the mass density, u> the frequency, the i 
component of the displacement fields u and C the fourth 
order stiffness tensor. This equation can be solved ana- 
lytically for a spherical system in the particular case of an 
isotropic stiffness tensor! 121] . In the case of an anisotropic 
stiffness tensor, we turn to a numerical solver. Using the 
scheme of Visscher et al.[29|, we develop the displacement 
field components on a polynomial basis. Solving Eq. (fljl 
then reduces to matrix algebra. We use in this study 
polynomials of order up to 20 to ensure the convergence 
of the method. Note that this parameter depends on the 
vibration mode of interest: high frequencies vibration 
modes commonly correspond to high spatial frequencies 
which require an extensive polynomial basis. 

2. Modes of interest and identification 

Severals experimental techniques enable the observa- 
tion of acoustic vibration modes, the main ones being 
inelastic neutron scattering(INS) 30], Raman scattering, 
time-resolved pump-probe experiments and far-infrared 
absorption. Whereas the absence of selection rules in 
INS makes it sensitive to all vibrations, the three other 
cited experimental techniques only probe a few acoustic 
vibration modes. Therefore, among the very numerous 
calculated vibration modes, we are mainly interested in 
those observable by Raman scattering experiments. Se- 
lection rules for Raman scattering and far-infrared ab- 
sorption for an isotropic solid sphere have been given by 
Duval[3l|. Only the breathing (spheroidal, I — 0) and 
quadrupolar (spheroidal, £ = 2) modes are Raman ac- 
tive. Note that time-resolved pump-probe experiments 
only detect the breathing modes of spherical NP [33] ■ 

The irreducible representations corresponding to all 
the vibration modes have been determined after the dis- 
placement fields have been calculated [15]. For a sphere 
made of a material having a cubic elasticity (such as silver 
or gold), the non-degenerate breathing mode transforms 
into an Ai g vibration (Oh point group). The quadrupo- 
lar mode splits into E g and T2 g vibrations [ill [13] and the 
corresponding frequency splitting is large making the us- 
age of elastic anisotropy a key factor for materials such 
as gold or silver. 



A. Elastic Calculation 

1. Solving the Navier equation with an anisotropic stiffness 

tensor 

The displacement fields u associated to each vibration 
mode are calculated in the scope of the linear elasticity 
theory by solving the Navier equation [28J: 



puj m + ^ Cijki a — q — = (1) 

jkl J 



B. Atomistic Calculation 

1. Calculation details 

We perform atomistic simulations using the EAM po- 
tentials developed by Clery and Roseto[33|. Table U re- 
ports structural and elastic properties calculated in the 
scope of this model for silver and gold using the program 
GULP P, HE [11. 

The calculated elastic coefficients have been used as 
input in the Navier equation (TTJ in Sect. Ill Al in order 
to focus on the comparison of the models and eliminate 
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silver gold 



a 4.078 A 4.079 A 

Cn 132.81 (131) 187.38 (187) 

C12 97.47 (97) 154.40 (155) 

C 4 4 51.11 (51) 44.71 (45) 



TABLE I: Cell parameter and elastic coefficients of Ag and An 
calculated using the EAM parameters of Cleri and Rosato[33l]. 
Elastic coefficients are expressed in GPa. Experimental 
data[13] for elastic coefficients are reported in parenthesis. 



possible differences arising from slightly different elastic 
tensors. 

Spherical NPs are designed with atoms initially placed 
on a perfect cubic close-packed crystal structure us- 
ing the bulk cell parameters of the modeled material. 
Our NPs do not specifically reproduced NPs with magic 
number(38j. In addition, our construction procedure in- 
evitably produces some steps and facets especially for 
very small NP. The total energy is then relaxed using 
a conjugate gradient algorithm: this procedure ensures 
that all vibration eigenfrcqucncies are real. Our model 
omits the possible existence of surface reconstructions 
and dangling bonds: these aspects could correctly be 
addressed using more sophisticated techniques (for in- 
stance DFT or tight binding). The dynamical matrix is 
then computed and diagonalized to obtain eigenfrequen- 
cies and eigenvectors. Eigenvectors are normalized for 
the usual scalar product in the real space of dimension 
3iV with N the number of atoms in the NP. In the fol- 
lowing, we study spherical NPs from 19 to 2072 atoms 
corresponding to diameter ranging from 1.6 nm to 4.3 
nm for silver NPs and from 10 to 2491 atoms correspond- 
ing to diameter ranging from 1.4 nm to 4.5 nm for gold 
NPs. The validity of the EAM potential for a 10 atoms 
NP is arguable and we admit that an ab-initio density 
functional calculation would be more appropriated in this 
case. However, in the following we compare the vibration 
properties of NPs on a large range of sizes. As a result, it 
does not seem judicious to use different atomistic mod- 
els for different sizes since it would also result in a slight 
modification of the elastic stiffness tensor. As a matter 
of fact, though quantitative results for the smallest NPs 
would need a more accurate modeling, we believe that re- 
sults described below are semi-quantitative in the sense 
that they give a good sketch of the vibration properties 
of NPs while reducing size. 

In the following, we will refer to atomistic modes for 
vibration modes calculated in the framework of the atom- 
istic EAM calculations. 



2. Modes identification 

While the linear elasticity theory assumes a continuous 
medium and thus provide a displacement field for each 
vibration mode, atomistic vibration modes are described 
by the data of each atom moves[24|]. For this reason, and 
because the precise symmetry of the NPs shapes is not 
fixed when relaxing the energy, the symmetry analysis we 
did in the case of the elastic calculation becomes more 
troublesome in the case of atomistic modes. Instead, we 
prefer to project elastic modes on atomistic ones as it has 
already been done in the case of semi-conductor NPs[23|. 
The scalar product of the pih elastic mode with the qth 
atomistic mode reads: 

~ -*elas \-*atomic ^ \ ^ ->elas ( r? \ s-?atomic /r>\ 

< U p \ U q >= 2-^ U P ( Ki > M q,i 

i£atoms 

where Ri are the atomic positions in the atomistic model 
and u a ^° mtc is the displacement of atoms i located at 

position Ri. Note that atomistic modes are normalized 
for this scalar product. We also normalize elastic modes 
so that: 

Y J <ul las \uf omlc > 2 =1 (3) 

A more relevant quantity is the sum of such squared 
projections of elastic modes onto atomistic ones for all de- 
generate elastic modes havin g th e same irreducible repre- 
sentation and eigenfrequencv jl5|| . More precisely, if £l e ias 
is a set of such elastic modes, we define its projection on 
atomistic ones as: 

/ \ 1 rTtelas I ^tatornic \2 

w? V N ^ 

where No elas is the number of modes in the set f2 e ; as . The 
denominator in Eq.((4]) ensures the normalisation condi- 
tion (similar to Eq.([3])) for the projection of a set of elas- 
tic modes. The quantity Q^" las does not depend on the 
choice of degenerate elastic displacements. As such it is 
suitable for the comparison with atomistic calculations 
and will be used in the following. 

III. RESULTS 

A. Elastic regime: Relevance of the anisotropy of 
the stiffness tensor 

1. Mode Projection 

When surface effects are negligible, i.e. the NPs size 
is larger than a characteristic size which will be estab- 
lished below, the elasticity theory reproduces very well 
the atomistic calculations. 
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Fig. [I] reports the projection of the set of fundamental 
quadrupolar E g and T 2g modes calculated by elasticity 
theory on atomistic modes for a 4.3 nm diameter silver 
spherical NP using Eq. (TJJ. In order to show the rele- 
vance of the anisotropy of the stiffness tensor in elastic 
calculations, Fig. [T] reports the projection of the set of 
fundamental quadrupolar elastic modes calculated using 
both an isotropic and an anisotropic stiffness tensor. The 
anisotropic elastic calculations use the stiffness tensor 
calculated from the EAM atomistic model and reported 
in Tab. HI The isotropic calculation uses a stiffness tensor 
obtained using the 3D averaged sound velocities. Atom- 
istic calculations intrinsically include the anisotropy of 
the elastic properties. 
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FIG. 1: (color online) Projection of the fundamental 
quadrupolar elastic vibration modes of a 4.3 nm diameter sil- 
ver spherical NP, calculated using the linear elasticity theory 
on atomistic ones using the scalar product given in Eq. ((2} 
and Eq. Q . Results are reported for both isotropic (above) 
and anisotropic (below) calculations. For isotropic calcula- 
tions, the quadrupolar elastic vibrations modes form a set 
of 5 degenerated modes with the same irreducible represen- 
tations. In the anisotropic calculations, quadrupolar modes 
splits in two sets of modes of respective symmetry E g and 
T2 g . Abscissa denotes atomistic modes using their frequency 
in cm -1 . 

In the isotropic elastic approximation the displace- 
ments of the 5 quadrupolar vibrations (m = 0, ±1,±2) 
correspond to a single frequency. In the anisotropic elas- 
tic approximation, this degeneracy is partially lifted. As 
a result, 2 of these "quadrupolar" vibrations have the 
same frequency (E g vibrations) and the 3 others have 
a different frequency (T2 g vibrations). It should be 
noted that the main effect of elastic anisotropy is to 
split the frequencies into 2 groups but that the corre- 



sponding displacements are almost unaffected p^j. Fig. Q] 
provides a simple way to check that taking into account 
anisotropy in the elastic calculations significantly im- 
proves the description. The quadrupolar modes from the 
isotropic elastic calculation have a significant projection 
onto atomistic modes over a much larger range than the 
corresponding anisotropic "quadrupolar" elastic modes. 
Therefore the anisotropic description is clearly an im- 
provement. The Eg vibrations have a significant pro- 
jection over a narrow frequency range. For the T 2g vi- 
brations, the projections are significant over a larger fre- 
quency range, but this is probably due to the light mixing 
which occur in the elastic calculations between the T 2g 
branches coming from the isotropic quadrupolar and tor- 
sional (I = 3) modes[15]. This mixing is different within 
the atomistic approach which explains the appearance of 
an additional peak at higher frequency (~ 16.3 cm -1 ). 

In the following, elastic calculations always use the 
anisotropic stiffness tensor defined from Tab. Q] 



B. Size dependence and mode number reduction 

Following the work of Cheng et al.[23| we report in 
Fig. [5] the evolution of the projection of the fundamen- 
tal breathing mode for different silver spherical NP sizes. 
This particular mode is the Ai g vibration with a fre- 
quency close to the isotropic breathing mode and hav- 
ing the largest volume variation [l5j. The projection of 
this elastic mode onto the atomistic ones is all the more 
peaked on a single atomic mode for big NPs. This obser- 
vation agrees with the expectation that elasticity works 
well for large NPs for which the continuous medium ap- 
proximation is more appropriated. Note that even for 
a very small NP, elasticity reproduces fairly well the 
atomistic modes. The same conclusion has been raised 
by Cheng et al.[23j while studying spherical germanium 
NPs. 

Reducing the size of the NP diminishes its number of 
atoms. Thus, it also reduces the number of vibration 
modes. From Fig. fJJ the breathing mode is evidenced 
even for the very small NP containing only 19 atoms. 
However, we expect that higher harmonics may disappear 
from the atomistic spectra while decreasing the NP size. 
In Fig. [31 we report the projection of the first overtone of 
the breathing mode (Ai g vibration with a frequency close 
to the first overtone of the isotropic breathing mode and 
having the largest volume variation) of silver spherical 
NPs on atomistic modes as a function of the NP size. 

Contrarily to the fundamental breathing mode, the 
projection of the first overtone spreads over several atom- 
istic modes for NPs made of 249 atoms and less. It is 
still possible to assign a main corresponding atomistic 
mode for a 249 atoms NP but no atomistic modes seem 
to describe the elastic mode for the 87 atoms NP. In- 
stead, this mode projects onto several atomistic modes 
with a weak scalar product. For smaller NPs, it mainly 
projects onto the same atomistic mode as the fundamen- 
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FIG. 2: (color online) Projection of the breathing vibration 
modes of silver spherical NP calculated using linear elastic- 
ity theory onto the atomistic modes using the scalar product 
given in Eq. <(2j) . The number of atoms in the NP is reported 
on each graph. The frequency in cm -1 of the atomistic modes 
is used on the abscissa. 
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FIG. 3: (color online) Projection of the first overtone of the 
breathing mode of silver spherical NPs calculated using linear 
elasticity theory onto the atomistic modes using the scalar 
product given in Eq. ([2j|. The number of atoms in the NP 



is reported on each graph. The frequency in 
atomistic modes is used on the abscissa. 
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C. Size dependence: breakdown of elasticity 
prediction 



1. Power Law 



tal breathing mode. This behavior is a manifestation of 
the disappearance of higher frequency elastic vibration 
modes due to the decrease of the number of atoms in 
the NPs. Indeed, while approximately 70 elastic modes 
have a frequency smaller than that of the fundamental 
breathing mode, this number increases to approximately 
500 for its first overtone. These numbers are to be com- 
pared with 3N — 3 which is the number of vibrations for 
a nanoparticle made of N atoms. We can crudely es- 
timate that the first overtone of the breathing mode is 
not well defined for N such that 3iV < 500 which is in 
agreement with the results presented in Fig. [3j The pro- 
jection onto the same atomistic mode as the fundamental 
breathing mode is due to the very close symmetry of the 
two modes. Note that in Fig. the abscissa scale extends 
over the whole spectrum of vibration states of the atom- 
istic model. Tamura and Ichinokawa [27| defined some 
maximum frequencies for spheroidal and torsional modes 
depending on the number of atoms in the NP in order 
to use the elastic model to calculate the specific heat of 
small NPs. A similar rule might apply to anisotropic cal- 
culations. But checking the rules proposed by Tamura 
and Ichnokawa from our calculations is beyond the scope 
of this paper. Current computer facilities make the calcu- 
lations of vibration modes of small NPs easy so that the 
rules proposed by Tamura and Ichonakawa is not crucial 
nowadays for interpreting experimental data. 



Using linear elasticity theory and from Eq. (Q]), the 
frequency of a given vibration mode of a solid sphere 
scales as 1/R where R is the radius of the sphere. This 
law which we will refer to as the Elastic Frequency Law 
(EFL) in the following, can be used to probe the validity 
of the elasticity theory. Recently, the validity of the EFL 
has been questioned for very small NP for which surface 
effects become non negligible (lil l25l|. 

In this section, we investigate the validity of the 
EFL in the case of metallic NP. We thus report the 
frequency of a given atomistic mode as a function of 
1/7V 1 / 3 (proportional to the radius of the NP) with N 
the number of atoms in the NP. This atomistic mode 
is defined as the one on which the projection of elastic 
modes is maximum ignoring the relative values of this 
projection compared to the other modes. Figure [4] 
reports the frequencies of four selected atomistic modes 
for spherical silver NP as a function of 1/iV 1 / 3 . We 
first choose four elastic modes: one E g mode and 
one T2 g mode coming from the fundamental isotropic 
quadrupolar mode and two Ai g modes coming from the 
fundamental and first overtone of the isotropic breathing 
mode. Then for each elastic mode and number of atoms 
N, we consider only the atomistic mode onto which 
the elastic mode has the largest projection and report 
its frequency as a function of 1/7V 1 / 3 . To identify the 
possible breakdown of the scaling law, we also plot the 
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FIG. 4: (color online) (above) Frequencies of the atomistic 
modes onto which the E g and T2 g modes coming from the 
fundamental quadrupolar modes and the Ai g modes com- 
ing from the fundamental (low frequency) and first overtone 
(high frequency) of the breathing mode preferentially project 
in spherical silver NPs as a function of the inverse of the num- 
ber of atom at the power one third. Dashed lines correspond 
to frequencies calculated by the elasticity theory: the num- 
ber of atoms is calculated from the ratio of the NP volume 
on the atomic volume computed from atomistic calculations, 
(below) Product of the frequency and the number of atoms at 
the power one third as a function of the inverse of the num- 
ber of atoms at the power one third. Linear elasticity predicts 
horizontal straight lines. 
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FIG. 5: (color online) (above) Frequencies of the atomistic 
modes onto which the E g and T2 g modes coming from the 
fundamental quadrupolar modes and the Ai g mode coming 
from the fundamental breathing mode preferentially project 
in spherical gold NPs as a function of the inverse of the num- 
ber of atom at the power one third. Dashed lines correspond 
to frequencies calculated by the elasticity theory: the num- 
ber of atoms is calculated from the ratio of the NP volume 
on the atomic volume computed from atomistic calculations, 
(below) Product of the frequency and the number of atoms at 
the power one third as a function of the inverse of the num- 
ber of atoms at the power one third. Linear elasticity predicts 
horizontal straight lines. 



product of the frequencies and iV 1 / 3 as a function of 
1/7V 1 / 3 . Note that the evaluation of the radius of a NP 
in the atomistic model may involve some technical issues 
regarding the precise definition of the surface, especially 
for small NPs: the surface is not exactly spherical and 
the position of the surface can be slightly different 
from the position of the surface atoms when taking into 
account atomistic radius. As a consequence, we prefer 
to report frequencies as a function of the unambiguous 
quantity 1/N 1 / 3 . 

Figure 0] shows that atomistic calculations agree well 
with the EFL for the atomistic quadrupolar modes and 
the atomistic fundamental breathing mode. Except for 



the smaller NP with 19 atoms, for which the frequency 
differs slightly from the EFL, all the other points follow 
the EFL. The frequencies of the atomistic fundamental 
breathing mode and the two fundamental quadrupolar 
modes differ by about 6% (8 cm -1 ) and 5 cm^ 1 and 
6 cm^ 1 respectively for 19 atoms NP. Difference for big- 
ger NP are smaller than 1 cm . 

Concerning the first overtone of the breathing mode, 
Fig. 2] shows that the atomistic frequency starts to differ 
significantly from the EFL for less than about 250 atoms 
corresponding to a characteristic diameter of 2.2 nm. 
This breakdown of the EFL could be related to surface 
effects flij (we exclude at this point the case of the 19 
atoms NP). Surface and surface relaxation effects thus 
do not affect equivalently the different vibration modes. 
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The relaxation of the surfaces spreads on a characteris- 
tic size A of a few Angstroms under the surface pH]. We 
can reasonably expect that modes with a small wave- 
length (such as harmonics of isotropic modes) are more 
affected than longer wavelength ones (or fundamental 
mode) by the presence of the surface. Obviously this 
claim also depends on the relative volume affected by 
the surface relaxation compared to the unaffected volume 
in the NP: frequencies of optical modes or border zone 
phonon modes in a semi-infinite medium are not affected 
by the presence of the surface. Finally, for the 19 atoms 
NP, the atomistic first overtone of the breathing mode 
is identical to the atomistic fundamental breathing mode 
because of the reduction of the number of modes and 
symmetry arguments as already discussed in Sect. IIII Bl 

In addition, we performed the same study for gold NP 
to investigate the dependence of our results on the mate- 
rial. Figure [5] reports the frequency of the fundamental 
quadrupolar and breathing atomistic modes for a spher- 
ical gold NP as a function of 1/A 1 / 3 . While the atom- 
istic quadrupolar modes reproduce the EFL fairly well, 
the atomistic fundamental breathing mode significantly 
shifts from the EFL at a characteristic size of about 87 
atoms corresponding to a diameter of about 1.6 nm: the 
frequency of the atomistic fundamental breathing mode 
differs by about 18% (21 cm" 1 ) from the EFL for 19 
atoms NP and by about 12% (8 cm -1 ) for the 87 atoms 
NP. As in the case of silver NP, we attribute the break- 
down of the EFL to surface relaxations effects. Note 
that our results suggest that surface effects seems to 
be more important in gold NPs than in silver NPs (for 
which the atomistic fundamental breathing mode repro- 
duce the EFL fairly well except for the smallest studied 
NP). This result is corroborated by the fact that surface 
relaxation in gold NP extends over a larger range than in 
silver [2ll [39j . and that the surface stress (to which sur- 
face relaxation is related) in gold is about two times the 
one in silver: this last result has been obtained for differ- 
ent surface orientations by atomistic simulations based 
on the EAM [II ELI and on the Modified EAM 0. 

In addition, to rule out the dependence of our results 
on the precise shape of the NP, we have also performed 
simulations on truncated cuboctahedron silver NPs: no 
significant deviation from the results presented here has 
been observed. 



IV. CONCLUSION 

By comparing atomistic and linear elasticity calcula- 
tions, we have exhibited the importance of the anisotropy 
of the stiffness tensor in elastic calculations. Especially, 
both atomistic and anisotropic elastic calculations lift the 
degeneracy of the quadrupolar modes. By projecting the 
elastic modes onto the atomistic ones, we evidenced the 
good agreement between the two models for the funda- 
mental breathing mode even for very small NPs. These 
two results suggest the applicability of the linear elas- 
tic theory to less than 2 nm diameter silver NPs. How- 
ever, we also evidenced the decreasing number of vibra- 
tion modes and the breakdown of the EFL when reducing 
the size. We conclude, that the applicability of the linear 
elasticity depends on the considered vibration mode and 
the NP size. Our results suggest that the breakdown of 
elasticity occurs for bigger NPs for high harmonics rather 
than for fundamentals modes. Moreover, critical sizes de- 
pend on the material. Such breakdown to our knowledge, 
has never been evidence experimentally in any material 
because of the difficulty to synthesize NPs of diameter of 
less than 2 nm with a narrow size distribution or to per- 
form single particle measurements on such objects. Our 
results suggest that the observation of the breakdown of 
the vibration frequency scaling in silver NPs will be haz- 
ardous, since it would imply some experiments on silver 
NP (or cluster) of about 20 atoms for the fundamental 
breathing and quadrupolar modes. This task though dif- 
ficult, seems more reachable in the case of gold NP since it 
would require measurements on NPs with less than about 
100 atoms. Finally, because the shape of very small NPs 
can hardly be described by a sphere, we have confirmed 
through simulations on truncated cuboctahedron silver 
NPs that our results do not strongly depend on the ex- 
act shape of the NPs. 
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